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Abstract 

OO ■ 

CN ' We show that the local magnetization in the massive boundary Ising model on 

■ the half-plane with boundary magnetic field satisfies second order linear differential 

43 . equation whose coefficients are expressed through Painleve function of the III kind. 

1 Introduction 

S: 

In the work [lj a very simple and elegant derivation of the famous Painleve equations for 
,— 1 . the spin-spin correlation function in the scaling Ising model with zero magnetic field was 
given. The approach used in that work was also applied in [2] to derive finite volume 
form factors of spin field in the Ising theory and in [3] to derive the differential equation 
for spin-spin correlation functions in the Ising theory on a pseudosphere. Here we ap- 
ply this approach to derive differential equation for local magnetization (i. e. one-point 
correlation function of spin field) in the boundary Ising model on the half-plane with 
boundary magnetic field. It turns out to be a second order linear differential equation 
whose coefficients are expressed through Painleve function of the III kind. Supplied with 
appropriate boundary condition it uniquely defines the local magnetization as a function 
^ of the distance to the boundary. Besides being interesting in itself, such a representation 

for local magnetization may be more convenient for numerical calculation in comparison 
with conventional form factor expansion, especially in the short distance region. 

As is well known [3J, there are two essentially different types of conformal boundary 
conditions (b. c.) in conformal Ising field theory. The so called "free" b. c. corresponds to 
the universality class represented by the lattice Ising model with unrestricted spins on the 
boundary. The so called "fixed" b. c. corresponds to the universality class represented by 
the lattice Ising model with boundary spins all fixed in the same direction ("+" or "— ", 
so there is more precisely two different "fixed" b. a). All this b. c. correspond to the 
fixed points of the boundary renormalization group flow. The most general local b. c. in 
the Ising field theory is the "free" b. c. perturbed by the boundary spin operator (which 
is the only non-trivial relevant boundary operator in the case of "free" b. a). This b. c. 
corresponds to the renormalization group flow from "free" b. c. towards one of the "fixed" 
b. c. More generally, one may consider conformal Ising field theory with "free" b. c, 
perturbed by both boundary spin operator and bulk thermal operator [5]. This theory 
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describes the continuum limit in the vicinity of the critical point of the lattice Ising model 
with zero magnetic field in the bulk and with boundary magnetic field being suitably 
rescaled. Let us briefly list known results about the local magnetisation in this theory 
defined on the half-plane. 

The form factor expansion for local magnetization a (t) was written down in [7j using 
exact expression for boundary state obtained in [6j: 



a 



(t) = (T exp ^XJ^/fcj 



(1) 
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where m ~ T — T c is the mass of a particle, h - scaling boundary magnetic field, y - 
distance from the boundary, cr - magnetization on the infinite plane, A = 1.28243 . . . 
is Glaisher's constant. Here and later on we always consider the low temperature phase 
T < T c , unless it is specially pointed out. It is also implied conformal normalization of 
spin field: 

i 

\x — x'\ 4 (a (x) a (x)} — > 1, as x — > x (4) 

With this normalization a (t, A) — > <jq as t — > oo. The expansion ([I])-© was first obtained 
in [9] from lattice model calculations. It was also shown in [8], [9] that in the cases of 
"free" (h = 0) and "fixed" (h — > ±oo) b. c. local magnetization can be expressed through 
Painleve function of the III kind: 



(Tfree (t) = CT exp <j \p it) + 1 
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cr fixed (t) = cr exp I -\(f (t) + ~ 



-^)-l + -(shV(r)-(^(r)) : 



l-e^ + -[sh 2 <f(r)-(<f' (r))" 



4 r w 4 

where ip (r) is the solution of radial sinh-Gordon equation: 



dr 



dr 



(5) 
(6) 



r 



sh 2ip 



satisfying asymptotic conditions: 
ip(r) = — In ( — -rQ 



0{r 4 Q 2 ) 



as r 



0, ft = In — r 



(7) 
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K (r) + (e _3r ) , as r -> oo 



(9) 



where 7 is the Euler's constant, i^o (^) is the modified Bessel function of zeroth order. 
As is known, (p (x) is related to Painleve function of the III kind 77 (x) as 77 (sc) = e"^^ 21 - 1 . 
More about this function see [TP], [IT]. 

In the case when bulk is critical (m = 0) it was shown in [I2j that: 

a (y) = h2-* ^ (2y)i * (l/2, 1, 87r/i 2 y) (10) 
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where 



oo 



y(a,c,x) = —— I e- xt t a ~ l {l + t) c ~ a - l dt (11) 

1 l (X ) 







is a solution of degenerate hypergeometric equation. 

The qualitative behavior of a (t, A) is well understood [?]. On the whole interval (0, oo) 
<?free (t) monotonically increases and af ixe d if) monotonically decreases, both approaching 
do as t — > oo. For small £0 

*/ree(*)~** (12) 
Ofixed{t) ~ r« (13) 

For < A < 2 a (t, A) remains monotonically increasing. Its values near the boundary are 
somewhat enhanced by the presence of boundary magnetic field, the leading term of its 
short distance asymptotic become dressed by logarithm: 

a(t)~t*kit (14) 

For A > 2 it possesses a maximum in some point. As A —>■ oo this maximum turns 
into a very sharp peak located in the region t ~ A -1 near the boundary, its shape being 
described by (TTOTh ffTTI) . For t < A -1 a (t, A) behaves as (fH}l . while for i > A -1 its behavior 
coincides with one under "fixed" b. c. (J6J. This dependence reflects the renormalization 
group cross-over between "free" and "fixed" b. c. 

The main result of this paper is that for arbitrary A: 

G (t, A) = U (t, A) a free (t) (15) 

where af ree (t) is given by (JBJ, and u (t, A) is the solution of differential equation: 

u" - (<f' - ch p + A) u' + ^A (if' - ch if + 1) u = (16) 
satisfying asymptotic condition: 

u (t) = 1+ o(t~5e-*V asi^oo (17) 



(Here (t) is the same function as in (jSJ), (jSJ) and the strokes stand for detivatives with 
respect to t.) 

Let us make some remarks on the equation (fl6|) . One can see that when A = 0, the 
only solution of (fT6|) satisfying (TlTj) is m (t) = 1. When A — >• oo (fTE|) turns into a first order 
differential equation which upon integrating and fixing integration constant with the help 
of ( TPTj) yields (EJ). In the massless limit (t — * 0, A — >• oo, £A kept fixed) (1T61) turns into a 
degenerate hypergeometric equation. Its solution can be fixed by "sewing" its asymptotic 
as tA — > oo with asymptotic of ([6]) as t — > 0, and this yields ffTUj) . ffTTl) . The fact that in 
massless limit we reproduce the result of [12] is not very surprising because the approach 
we used to derive ffTBI) . ([TBI) is a generalization of one used in [12]. Concerning the relation 
between the form factor expansion ([T]), ([2]), ([3]) and our result ffTBI) . ffTBI) we just note that 
it seems to be very difficult to show directly that f|l]), f|2]) satisfy ffT5|) . ffT6|) . In any case, 
it is beyond the analytic abilities of the author. 



^^Note that comparing the coefficient in the short distance asymptotic of <Jfi xe d{i), that follows 

from |6]), with the result <Tfi X ed(y) = 2^ (2y)~ s (for m = 0) of 13J , one obtains the identity 
J oo (l_ e - V W) dr . = l n2 
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Being second order linear differential equation, (|T6|) possesses two linearly independent 
solutions. Their asymptotics as t — > oo are u\ (t) ~ 1 and 112 (t) ~ e*- Hence, for 
A > 1 the condition that u (t) — > 1 as t — > 00 is sufficient to fix the solution uniquely. For 
A < 1 more strict condition (fTTj) is required, which follows from form factor expansion 
([1]), ([2]). Another linearly independent solution in this case also has physical meaning. 
As explained in [6], for A < 1 there exists metastable state characterized by asymptotic 
behavior a (t) — > — <jq as t —>■ 00 and corresponding to the boundary bound state in the 
hamiltonian picture with "space" being half-line and "time" axis being parallel to the 
boundary. The local magnetization o\ it, A) in this state can be obtained from a (t, A) by 
analytic continuation h — > — h. Clearly, it is also a solution of (fT6j) . As it was shown in 
[H] its asymptotic as t — > 00 is: 

1 

u x (t, A) = -a + a (^a) ' e~ (1 - A)t + ^= - l) ri«T* + (r§ e"*) (18) 

The presence of exponential term ~ e~ ( - 1 ~ A - ) * in (Tl8|) agrees with ([TBI) . 

In the rest of the paper we present the details of our derivation of (TT5T) . (TT61) . 



2 Ising field theory in the bulk 

In this section we briefly recall some well known facts p] about the structure of the Ising 
field theory in the bulk needed for further computations. As is known the Ising field theory 
in zero magnetic field is equivalent to the free Majorana fermion theory with euclidean 
action: 

S = — (tpdtp + rjjdrjj — irmpip) d 2 x (19) 



2tt „ 

Here we have assumed that the theory is defined on an infinite plane M 2 , whose points x 
are labelled by cartesian coordinates (x, y) = (x (x) , y(x)), and d 2 x = dxdy. Complex 
coordinates are defined as z (x) =x + iy, z (x) =x -My, and the derivatives d,d in f fT9"j) 
stand for d z = \ (<9 X — id y ) and d% = \ (<9 X + id Y ) respectively. The ciral components if), if) 
of fermi field satisfy Dirac's equations: 

= di) = ^ (20) 

Their normalization in the action ( fT9|) corresponds to the following short- distance limit 
of the operator products 

zif>(x)il>(p) zi>(x)i>(p)->l, ass^O (21) 

The order o (x) and disorder fi (x) fields are semi-local with respect to the fermi fields; 
the products 

xf)(x)a(0), ip(x)fi(0), ]>(x)<t(0), ]>(x)h(0) (22) 

acquire a minus sign when the point x is taken around zero point. The fields if) (x) and 
if> (x) in the products (1221) can be expanded in the complete set of solutions of Dirac's 
equations ( 1201) having this monodromy property: 



4 



where 



S „W;-V2J r i"+2) -ie'Ki)»/":, M 1 <24) 



f n Ix) \ fm\2 



i 2 

1" 



v 2 y r T + 2 J e -*MW , fwl 1 (25) 



n "2 



(here r, # are polar coordinates, i. e. z — re te , z = re~ td and I v are modified Bessel 
functions). The coefficients a n , a n in (1251) are understood as operators acting on the space 
of fields. 

It can be easily shown that for any two solutions = Uj) x (x) , ^ (x)), ^2 = 
(■^2 j 1 02 °f Dhac's equations which change sign after the point x is taken around 
zero the integral 

* 2 ) = 77- / V'i 0*0 ^2 0*0 ^ - $1 ( x ) ^2 0*0 ^ (26) 



2m 

Co 

over a contour Cq encircling zero (in counter-clockwise direction) does not change under 
continuous deformation of Cq and therefore defines a bilinear form on the space of such so- 
lutions. The solutions U n = (u n ,u n ) and V n = (v n ,v n ) satisfy the following orthogonality 
properties with respect to this bilinear form: 

(U n ,V m ) = (27) 

Let us also write down the following differentiation formulas which are useful in compu- 
tations with U n and V n : 

( 1\ - m 2 

dU n = ^ n --jU n ^, dU n = 2{2n + i) U n+1 

m 2 ( l\ 

dv " = 2|aTTT) l/ " +1 ' dv " = (" " 2) v -' (28) 

(here we denote dU n = (du n (x) , du n (x)), etc.). 

Using relations (T271) one can express operators a n , a n in terms of contour integrals: 

a n = — : ® u n [x) i[> (x) dz — u n (x) (x) dz (29) 



2m 

c 

a n = j) v n (x) ip (x) dz - v n (x) i) (x) dz (30) 

Co 

This representation can be used to show that they satisfy canonical commutation relations: 

{a n ,a m } = (31) 

The fields a and /z are "primary" with respect to the algebra ( 13T1) . i. e. they satisfy 
relations: 

a n cr = , a n a = , a n jj = , a n jj, = (32) 
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for n > 0, as well as 



UJ UJ 



a a = -j=n , a /i = -j=o 

_ UJ 

a cr = — n , a fi = —a (33) 

where uj = e i7r//4 and uj = e~ m l 4 . These equations define the fields a and /i up to normal- 
ization. In what follows we will assume conformal normalization of fields a and \i: 

i i 
|x| 4 (7 (x) (J (0) — > 1 , \x\ 4 fi (x) fi (0) — > 1 , as x — > (34) 

As it is shown in [1], first and second order descendants of a and [i with respect to the 
algebra a n , a n are expressed in terms of coordinate derivatives of a and fi: 

CL-\0 

a_ia = —=4:dfi , a_i// = —=4:da (35) 







= 


^4da 
V2 




cL 


-i/i = 


V2 


UJ 8 2 
v/23 ^' 


a_ 


-2A* = 


Q 8 2 
—=-o a 
v/23 


w 8^ 
v/23 ^' 


a_ 


-2/i = 


^8-2 

—=-o a 

v/23 



(36) 

This observation is very important for the method of [1] to work. 

The Majorana theory ( fl9l) corresponds to both high and low temperature phases of 
the Ising model in the vicinity of its critical point T c depending of the choice of the sign 
of the mass parameter m in ( fl9i) . Our definition in ( 1331) corresponds to the identification 
of the case m > with the ordered phase T < T c , while the case m < is identified with 
the disordered phase T > T c . From now on we will consider the ordered phase m > 0. 



3 "Free" and "fixed" boundary conditions 

In this section we consider the Ising field theory defined on the half-plane y > 0. As a 
warm-up exercise let us first rederive formulas ([5l) ? for local magnetization in the cases 
of "free" and "fixed" b. c. Explicit expressions for boundary states for "free" and "fixed" 
b. c. was obtained in [3]. It follows from this expressions that the fields ip and ip satisfy 
the following b. c: 

(V>--0)| y=o = (for "free" b. c.) (37) 
(ip+^)\ y=Q = (for "fixed" b. c.) (38) 

Suppose that (x (x) ,x( x )) i s a double- valued solution of Dirac's equations defined on the 
half-plane y > (with punctured point xq) such that it changes sign when x is taken 
around Xo, decays sufficiently fast as \x\ — > oo, and satisfies b. c: 

(X-X)\ y=0 = (for "free" b. c.) (39) 
(X + X)l v=0 = (for "fixed" b. c.) (40) 
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Then the following identity holds: 



X (x) if) (x) dz — x (x) if) (x) dz ) \i (xq) ) = 



(41) 



where C Xo is a contour encircling the point xq. This is because the integral on the left- 
hand side of (JUJ) does not changes under continuous deformations of the contour C XQ and 
therefore one can deform it in such a way that it constitutes of two parts: which 
tends to infinity and C& which passes along the boundary. Then the integral along 
is zero because % ( x ) an d X ( x ) decay at infinity and the integral along Cb is zero due to 
(1391) or ( l40l) . On the other hand one can shrink C xo to a small circle around the point 
Xq and express the left-hand side of f|4Tl) in terms of descendents of \i using the operator 
product expansions of if>(x), if) (x) with \i (x ). If the singularity of (% (x) ,x( x )) a ^ the 
point Xq is not too strong only descendents of not higher than the second order appear 
in this expression. Since they are expressed in terms of coordinate derivatives of a this 
will lead to a differential equation for local magnetization (a(x)). The only problem is 
to find a solution of Dirac's equations (x (x) , x (x)) satisfying the above conditions. For 
this purpose we will use the following trick. We will search for this solution in the form 
of the linear combination: 



x( x ) 



(ip(x)a{x )n(Px )) 
' 1 (ip(x)a(x Q ) fi(Px )) Q 



+ c 2 d 



+ cd( ^( x ) a ( x o) »(Pxo)) Q 
3 ^ (ip (x)a (x ) fi(Px )) 



+ c 4 



{ip (x) a (x ) fx (Px )) 
(if; (x)(7(x ) fi(Px )) Q 

(ip (x) /I (x ) a {Px )) 
(if>{x) /i(x )a(Px )) 



+ 



+ 



+ c 5 <9 



(ip (x) n(x )a(Px )) 
(ip(x) /j(x )cr(Pxo)) 



+ cd( (f( x )^( x o) a ( Px o)) 
6 y (ip(x) /j(x )a(Px )) 



(42) 



where (. . .) denotes a correlation function in the Ising field theory defined on the plane 
and P denotes the reflection in the line y = (i. e. P (x, y) = (x, — y)). Obviously 
each term in (142p is non-zero and function of X IS cL solution of Dirac's equations, 
change sign when x is taken around Xo and decay at infinity. The coefficients in this linear 
combination can be determined from the requirement that it satisfies (|39|) or (140]) . Note 
that we do not need to know the functions x ( x ) an d X i x ) i n explicit form. What we 
really need is several terms of their short- distance asymptotics as x — > xq, but the latter 
can be expressed in terms of the two-point functions (a (xq) a (Pxq)) = G (2my ) and 
(fi (xq) fi(Pxo)) = G(2my ). As is well known \]Xj\ (see also pQ) this functions can be 
expressed in terms of Painleve function of the III kind as follows: 



G (t) = (T ch ( -ip (t) ) exp 



1 f°° ( 2\ 

-J r ysh 2 ip (r) — (ip' (r)) j dr 



(43) 



G(t) = a Q sh 



exp 



r ( sh 2 ip (r) — (if' (r)) 2 ) dr 



where ip (t) is the same function as in (j3J), (jSJ). 

Under parity transformation P fermi fields transform as: 



(44) 



P 



—iif> 
iif) 



(45) 
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We have therefore the following identities which follows from the invariance of correlation 
functions under parity transformation: 

(ip{x)a{x ) p(Px )) \ y=0 = i(ip{x) n{xo)(r(Px )) \ y=0 (46) 

(ip(x)fi(x )a{Px )) \ y=Q = i(ip(x)a(x ) fi(Px )) \ y=Q (47) 

(8$ (x) a (x ) //(Px )) | y=0 = i(dip(x) p(x ) a (Px )) \ y=0 (48) 

(d4)(x)n(x )a(Px )) \ y=0 = i(dip(x)a(x )fi(Px )) Q \ y=0 (49) 

It follows from this identities and Dirac's equations that from all correlation functions that 
present in the expression (1421) only four are linearly independent functions of x on the line 
y = (for example (ip (x) a (x ) [i (Px )), (ip (x) /J, (x ) a (Px )) , (dip (x) a (x ) {Px )), 
and (dip (x) fi (x ) a (Px ))). Hence requiring that (j4~2l) satisfy fl39|) or fj40|) one obtains 
four linear constraints for six coefficients ci, . . .c§. It turns out that they have non-zero 
solutions. One of the solutions corresponds to the following linear combination (it does 
not matter what of the solutions to choose): 

(ip (x) a (x ) \i {Px )) \_ w ( (ip (x) a (x ) fj, (Px )) \ 
{ip (x) a (a? ) H {Pxa)) J \ (V> ( x ) a ( x o) V ( p ^o)) J 

mf (ip(x)/i(x )a(Px )) \ (ip (x) /i (x ) a (Px )} \ . , 

+ 2 \& (x) /i (x ) a (Px )) Q J a \(^(x)n (x ) a (Px )) J ' [M) 

for "free" b. c, and 

X(x)\_.m/ (iP(x)cr(x )fi(Px )) \ _ d f (ip(x)a(x )n(Px )) . 



x 


(x) 


\ ml 


V x 


(*) . 


I = ~2 \ 



/ 2 V (^(a:)o-(xo)/i(Pxo)) ^ \ (^( I ) ff ( a: o)/i(^o)) 
m/ (x) // (x Q ) (7 (Px )) \ (^(z)/i(xo)«7(Par )) 



+ 2 { (iP (x) p (x ) a (Px ))° J ^ V <^ (*) ^ (*o) ff ( Px ^i ) ' (51) 
for "fixed" b. c. 

It is now straightforward but somewhat tedious exercise to substitute (|50|) and (|5"T|) 
in (14 ip and evaluate the left-hand side. One has to expand ip and ip using (1231) . than to 
evaluate contour integrals using (1271) . (f2"8"l) and than to act by the operators a n , a n on a 
and using (1321) . (1331) and (135|) . Due to (132]) all terms with descendants of higher than 
the first order vanish. Finally, taking into account that (a (xq)) = a (2my ) depends only 
on yo due to translation invariance, one obtains the following differential equations: 



2 (G - Gj a' free - [G 1 -G' + Gj a free = (52) 

2 (G + G) a' flxed - (g 1 + & + flxed = (53) 

(the stroke denotes derivative with respect to t — 2my ). Integrating this equations, sub- 
stituting (|43J) . (jUj) and fixing integration constants with the help of asymptotic condition 
er (t) — > cr as t — ► oo one obtains (jSJ) and (jHJ). 

Let us now consider the high temperature phase T > T c . The differential equations in 
this case can be obtained from (152"j) . (1531) by substitution m — » — m, G ^ G: 

2 (G - G) a' /ree - (G' - G' + G) a free = (54) 
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2 (G + G) cr^ ed - (C + G" - G) <7 /feced = (55) 

The only solution of ( 1541) that does not grow exponentially as t — > oo is a f Tee = 0, while 
from ( 1551) we obtain: 

& fixed, T>T C = £ 2t (T fixed, T<T C (56) 

in agreement with [9]. This confirms our identification of the case m > with the low 
temperature phase. Had we chosen the other choice, we would obtain the exponentially 
growing solution for a fi xe d in the high temperature phase. 



4 Boundary magnetic field 

Let us now consider the general case of "free" b. c. perturbed by boundary spin operator 
ob- The latter is identified with degenerate primary boundary field with dimension A = 
1/2 [4 J . It can be written in terms of fermion fields as follows [6] : 

o B (x) = io (x) (ip (x) + # (ar)) | y=Q (57) 
where a (x) is additional fermionic degree of freedom with two-point function 

(a(x)a(x')) /ree = ^sign(x-x') (58) 
The action of the theory has therefore the following form: 

oo oo 

S = — I dx I dj (?pd?p + ^)d^) — im^ipj + 



oo oo 

+ j (^-^ (^)\y=o + l ad - a ^ dx + ih j (ij + ij)\ y=0 dx (59) 

— oo — oo 

It leads to the following b. c. for fermion fields [6]: 

^&-ty\ y=0 =-im\^ + ty\ y=0 (60) 

where A = 4irh 2 /m. We can now proceed in the same way as in the previous section 
but now instead of f[59"j) or fT4Dl we should require the functions x an( i X to satisfy the 
condition: 

d 

[X - X) l y=0 = ~ imX [X + X) | y=0 ( 61 ) 
in order to write down the Ward identity (|4ip . It turns out that in this case it is necessary 



to include also terms with second order derivatives in the linear combination 042 1) in order 



9 



to satisfy flSlj) . As a result one obtains the following linear combination: 



2A) 



(tp (x) a (xq) /i (Pxq)) 
tp(x)a(x ) n(Px )) Q 



m 



2 

+ td 2 



I — O 

2 



l-2\)d 



(tp (x) a (x ) [i {Px )) 
(tp (x)a(x ) fi(Px )) 



m 5 

a 

2 



(tp (x) a (a? ) n (Px Q )) 
(ip(x)a(xo) fJ l (Px )) 

{ip (x) /J,(x )cr(Px )) 

(ip(x) /i(x )a(Px )) 



rn\ 2 , 
(y) (1-2A) 



777 — 

•j(l-2A)0 



+ 5 2 



(tp (x)a(xo)fi(Px )) Q 
(ip(x)a(x ) n{Px )) 

(tp (x) n(xo)a(Px )) 
(tp (x) n(x )a(Px )) Q 

(tp (x) fx (x ) a (Px )) 
(tp(x) ii(x )cr(Px )) 

(tp(x) fi(x )a(Px )) 
(tp(x) fi(x )cr(Px )) 



+ 



+ 



(62) 



Substituting it in (f4T!) and evaluating the left-hand side we obtain the following differential 
equation for local magnetization a (t) : 



G + G) a" 



G' + G'-G + XG + G 



1 

+ 4 



G" + G" - - ( G' + G' ) - 2G' - G + 2A ( G' + G' + G 



a = (63) 



It can be brought to a simpler form ffl6|) by means of substitution ffl5|) . 
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